Perfect state transfer in XX chains induced by boundary magnetic fields 
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A recent numerical study of short chains found near-perfect quantum state transfer between the 
boundary sites of a spin- 1/2 XX chain if a sufficiently strong magnetic field acts on these sites. We 
show that the phenomenon is based on a pair of states strongly localized at the boundaries of the 
system and provide a simple quantitative analytical explanation. 



PACS numbers: 03.67.Hk,03.65.Yz,75.10.Pq,75.40.Gb 



I. INTRODUCTION 

The transfer of quantum states between different ele- 
ments of a quantum computer is an important and chal- 
lenging task in quantum information processing. Chains 
of interacting spin- 1/2 particles were proposed as a means 
to achieve this task. Many different types of quantum 
spin chains have been discussed, as evident, for example, 
from the reviews [l| and[3- In the most frequently studied 
scenario the first site of a finite spin chain is prepared in 
a definite state which is then supposed to be transferred 
to the last site by the natural Hamiltonian dynamics. 

Since the elementary excitations of spin chains are of- 
ten dispersive, the initially sharply localized quantum 
state broadens as it propagates and thus quantum infor- 
mation degrades^ so that the state received is not iden- 
tical to the state transmitted. In view of limits to fault 
tolerance in quantum information processing it is, how- 
ever, desirable to achieve perfect state transfer (PST). 
Engineered chains were shown to be attractive means in 
reaching this go&}^^—. In these chains, all spin-spin cou- 
plings and / or local fields are set at certain values which 
guarantee PST. 

However, from a practical point of view, it would be 
more desirable to achieve PST by simpler means, for ex- 
ample by manipulating only 0(1) system parameters in- 
stead of the 0{N) coupling constants which must be ad- 
justed for an TV-site engineered PST chain. This line of 
thought has been pursued in a number of studies in which 
the ends of a homogeneous spin chain, i.e., one with con- 
stant couplings, were connected to a single transmitting 
and receiving qubit, respectively, by adjustable coupling 
constantsiSrJ^ which infiuence the transfer of states via 
different mechanisms, depending on other details of the 
systems employed. 

In a similar vein, Casaccino et alJ^ suggested to use 
sufficiently strong magnetic fields acting exclusively on 
the boundary spins of an otherwise homogeneous spin 
chain. In their numerical results, they indeed observed 
near-perfect state transfer, however, with the transfer 
time growing strongly with both the chain length and 
the magnetic field strength. 

In this paper, we explain the effect observed in Ref. 
[l6l by analytical considerations. We show that bound- 
ary magnetic fields of sufficient strength h produce ex- 



ceptional energy eigenstates of the chain. These states 
are exponentially localized at the boundaries with local- 
ization length ^ (ln/i)~^. The PST is completely de- 
termined by these states and the dynamically relevant 
energy differences become exponentially small as the sys- 
tem size grows. Correspondingly the time scales become 
exponentially large, making this kind of PST unfortu- 
nately impractical for long chains. Furthermore, since 
the PST depends on an exponentially weak link between 
transmitter and receiver (the central region of the chain, 
where the amplitudes of the relevant states arc exponen- 
tially small) one might speculate that any slight imper- 
fection is bound to strongly affect the PST performance. 
The robustness of state transfer against disorder in the 
couplings was recently studied for several kinds of both 
engineered and boundary-controlled spin chain a^^'^^ . 

In Sec. In] of this paper we present analytical solutions 
for the energy eigenstates and eigenvalues of the homo- 
geneous semi-infinite chain with a magnetic field at the 
first site. The results thus obtained suggest a simple per- 
turbative calculation and an approximate expression for 
the exact result for the finite chain with boundary fields 
at both ends, which we discuss in Sec. IIII | a nd which 
explains the numerical observations in Ref. Ila . 



II. THE SEMI-INFINITE CHAIN: 
RESULTS 
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The spin- 1/2 XX chain with boundary magnetic fields 
is defined by the Hamiltonian 
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where S'f , S'f , and S'f are the usual spin- 1/2 operators 
with eigenvalues ±1/2. The most popular scenario in 
quantum information transport studies uses an initial 
product state, = 0)) = \tp) ®f^2 I with \tp) = 
a I t) i)i i-e-i a superposition of the state with no spin 
up and a state with one single spin up located at site 1. 
Since the Hamiltonian conserves the number of spins 'up' 
and since the zero spin up component does not evolve in 
time it suffices to follow the evolution of the single spin 
up component. The corresponding A^-dimensional sub- 
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space is spanned by the states \i) = S^\ || ... |) and the 
Hamiltonian matrix Hij is given by 
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For h = the eigenvectors u'^ of the matrix ([2|) are given 
by their components 



N + 1 
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with 1/ = 1, A^, and the energy eigenvalues are 

VK 



E[y = 2 J cos 
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By the Hellmann-Feynman theoremi^r— all eigenvalues 
grow if a magnetic field /i > is switched on. To see 
more precisely what happens for h ^ we now consider 
the case N ^ oo, i.e., we effectively study a semi- infinite 
chain with boundary magnetic field h. Since we are look- 
ing for states localized near the boundary / = 1 we use 
the Ansatz 

ui = e"''' (5) 
which turns out to be an eigenvector for h > J, with 

K = In (6) 

and 



Obviously 6 has a singularity for /i > J at some q. How- 
ever, that singularity is only a jump from tt/2 to — 7r/2 
(or vice versa) corresponding to an irrelevant sign change 
in the eigenvector components. 

If N is finite (but large) the amplitudes of all extended 
eigenstates scale as N~-^/^, while the amplitude ^ of 
the localized state at the boundary site is constant and 
approaches unity as h grows. This is the key to the 
(quasi-) PST reported in Ref. [l^. As N or h grow the 
localized states at the boundaries will increasingly dom- 
inate the transfer of localized quantum information be- 
tween the ends of the chain, and we expect the speed 
of that transfer to be proportional to the overlap be- 
tween the states localized at the left and right ends of 
the chain, respectively, which by ([S]) and ^ is propor- 



tional to e 



III. THE FINITE CHAIN: APPROXIMATE 
RESULTS AND PERFECT STATE TRANSFER 

We now consider the full model ([T|), with a strong 
boundary magnetic field ft, ^ J at both ends of the A^- 
site chain. We then have two extremely localized energy 
eigenstates with energies far above all other states. The 
system can thus be approximately described by two iso- 
lated spins in a strong field and an {N — 2)-site chain 
in between. The corresponding Hamiltonian matrix is 
obtained from ^ by setting 



E = 2 J cosh K = h + —. 

h 



(7) 



We thus find a critical field value, above which one energy 
eigenstate splits off from the quasi-continuum (j4|) and 
moves upward with growing h. The localization length 
of that state is proportional to {Inh/ J)~^ and its ampli- 
tude at the boundary taking into account normalization 
is given by 



("i)^ = l-(^ 



(8) 



The remaining eigenstates of the system are extended 
and similar to ([3]); the appropriate Ansatz is 



ui = sin{ql + S), 



(9) 



where g is a wave vector between and tt. For finite N 
the value of q is fixed by the boundary condition at the 
far end of the system, leading to g = vir/ (N+l) for ft = 
as in Eq. (|3]), while for iV — > oo g becomes continuous 
and can take any value between and tt. The state ([9]) 
yields the eigenvalue 



E = 2 J cos g 
and the phase shift 5 obeys 

ft sin q 



tan (5 



J — h cos q 



(10) 
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Hi2 — H21 — Hn^n^i — H 
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= 0. 
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Taking parity into account, the two localized states are 
given by the amplitudes 



ut = = = otherwise (13) 
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both states have energy E± = ft. The N ~ 2 extended 
states are given by 
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0; < 
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(i/ = l,...,A^ — 2), with energy eigenvalues 



Eu = 2 J cos ■ 



N - 1 



(15) 



Returning to the full problem defined by the Hamil- 
tonian matrix ([2]) we expect modified energy eigenvalues 
E± for the formerly degenerate isolated boundary states 
|±). The initial and final states of the desired perfect 
state transfer can be written in terms of the states ([T^: 



|l) = ^(l+) + |-)) 



1 



l^) = ^(l+)-|-)) (16) 



and consequently the transfer time r is equal to the time 
during which a phase difference tt develops between the 
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energy eigenstatcs |+) and |— ), that is 



(17) 



The energy difference — can be determined ap- 
proximately by treating the matrix elements J between 
sites 1 and 2 and sites N — 1 and N, respectively, as 
perturbations. The corresponding perturbation opera- 
tor V then is a X iV matrix with four nonzero entries 
J. The matrix elements of V between the unperturbed 
eigenstates vanish and thus first-order degenerate 
perturbation theory does not lift the degeneracy. In sec- 
ond order, the corrections to the energy are given^^ by 
the eigenvalues of the matrix 



^ {n\VW){i.\V\n') 



E, 
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where |n) and \n') are the unperturbed degenerate energy 

eigenstates with eigenvalue En'' = h, and {v) runs over 
all other unperturbed eigenstates with eigenvalues 

Ei°^ (HID. In our case the degenerate subspace is spanned 
by 1+) and |— ), or equivalently, by |1) and \N), compare 
(jl6p . The matrix elements of the perturbation in the 



basis {|1), I A'')} are (compare ((T4)) ') 
{1\V\,^) = Ju^, = J ^ ^ 



N - 1 



N - 1 



(19) 



{N\VW) = Ju-^^,, = {-ir+'{l\VW). (20) 



and 



The diagonal elements {n = n') of the 2x2 matrix ([T 
thus turn out to be equal, leading to a mere energy shift 
The non-diagonal elements are both equal to 



N-l 



{l\V\v){v\V\N) _ 2P (-1)^+1 sin^ 

1. /,_2Jcos^ - iV-1 /i-2Jcos^ ' 

(21) 

leading to an energy splitting equal to twice the non- 
diagonal element (PT|) : 



E+~E^ = 



N 



AJ {^^ (-l)-+^sin^^ 
1^1 7-2cos^ 
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An asymptotic expansion of (j22p in the small parameter J /h reads 



E. 



E^ = 2 J 
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This confirms the behavior which was to be expected 
based on the results of Sec. [Tl] and it explains the expo- 
nential growth of the transfer time p7|) with both N and 
h/J. 

The exact expression of the energy difference Ej^ — E^ 
can be found and approximated by the following sim- 
ple consideration. The boundary-localized symmetric 
and antisymmetric eigenstates are determined using an 
ansatz of the form 

uf = e-^'±e-<''+^-'\ (24) 

The corresponding energy eigenvalues are again given by 

S± = 2JcoshK±, (25) 

(compare Eq.©) but now the values k± fulfill the tran- 
scendental equation 

e^^i = ± e-^«± (^^e«± - 1^ (26) 

which results from the condition for the eigenvalue at 
sites 1 or A''. This equation can be solved itcratively. 



starting from kq = In -j (compare Eq.®). To lowest 
order in {J/h)^^"^ this yields 

^-^-"(f) ('-(() )+°((S )■ 

For comparison to the numerical results of Ref. [161 we 
identify h with A^; from that reference and set J = 2. 



h 


Ref. 16 


exact diag. 


perturbative 


order (J/h)^-^ 


10 


99 


107 


90 


107 


20 


8010 


801 


770 


801 


30 


2665 


2674 


2627 


2674 


40 


6260 


6315 


6252 


6315 


50 


12294 


12311 


12233 


12311 



TABLE I: Transfer times for a chain of length N = 5 with 
different values of the boundary field h as obtained by various 
methods. See text for details. 

In Table [J we compare the transfer time for N ~ 5 as 
obtained by different methods and for different values of 
h. The column labelled "Ref. contains the numbers 
reported in that reference (with an obvious typing error 
for h = 20). The next column shows the transfer time 
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(fT7|) from the difference of the two dominant eigenvalues 
£'+ and as obtained from an exact numerical diag- 
onalization of the Hamiltonian ([2]), whereas the fourth 
column is obtained from the perturbation calculation de- 
scribed above. 

The asymptotic expansion ([23]) agrees to the expres- 
sion (|22p to a much higher precision than that of the num- 
bers in the table. However, it can be seen that the per- 
turbative result approaches the exact one only for fairly 
large values of h. On the other hand, the lowest order in 

(^) , given in Eq. (P7)) . agrees with the exact values 
to all digits given (compare the last column in Table H]). 

We have performed similar comparisons for chain 
lengths N = 7 and = 12, obtaining similar results, 
that is, for sufficiently strong boundary fields the pertur- 
bation calculation yields the transfer time to a precision 
of one percent or better while the lowest order fits the 
exact result within 10~^ for iV = 7 and within 10~^^ for 
N = 12, as was to be expected from Eq. ^7} . However, 
at = 12 the difference between the two leading en- 
ergy eigenvalues approaches the precision limits of stan- 
dard numerical procedures and the corresponding trans- 
fer times become prohibitively long by any standards. 

Finally, we briefly discuss the degree of fidelity which 
can be reached by the scheme discussed here. A reason- 
able measure of the transfer fidelity between the states 
|1) and |iV) is 

hN{t) = \{l\e~^"'\N)\\ (28) 

Employing the energy eigcnstates | v) and the eigenvalues 
we can write 

(l|e-^*|7V) = ^(l|^)(HiV)e-^-* (29) 

= (l|+)(+|Ar)e-*^+* 
+ (l|-)(-|iV)e-^-* 
+ ^(l|^.)(^|iV)e-^-*, 



where we have separated the dominant (localized) eigcn- 
states |±) from the remainig (extended) eigenstates 
each of which yields a contribution 0{1/N), with a quasi- 
random phase. The transfer time r is the instant when 
the two dominant terms in (j29p interfere constructively, 
maximizing fiN and thus leading to (compare ([8])) 

|(l|e-*^n^)l«^l- (30) 

or 

/m(r) = 1-2(^0'. (31) 

The dominant time dependence of fiN{t) is a slow os- 
cillation between zero and the maximum value pi|) with 
period 2t, with superimposed rapid small-amplitude os- 
cillations from the subdominant terms in (|29p . This is 
precisely what numerical calculations for small N show. 

The value /ijv (t) does not depend on the chain length 
N and it grows (increasingly slowly) as h grows. Note, 
however, that the transfer time r grows exponentially as 
h grows. For a given N the choice of h thus is the result 
of a trade-off between transfer time ([T7|) and fidelity ([3T|) . 
depending on the relative importance assigned to each of 
these two figures of merit. Hence, our analytical results 
for transfer time and fidelity provide the basis for finding 
an optimum trade-off in practice. 
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